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In recent years, there are attempts to describe quantization of planetary distance 
based on time-independent gravitational Schrédinger equation, including Rubcic & 
Rubcic’s method and also Nottale’s Scale Relativity method. Nonetheless, there is 
no solution yet for time-dependent gravitational Schrédinger equation (TDGSE). In 
the present paper, a numerical solution of time-dependent gravitational Schrédinger 
equation is presented, apparently for the first time. This numerical solution leads 
to gravitational Bohr-radius, as expected. In the subsequent section, we also discuss 
plausible extension of this gravitational Schrödinger equation to include the effect 
of phion condensate via Gross-Pitaevskii equation, as described recently by Moffat. 
Alternatively one can consider this condensate from the viewpoint of Bogoliubov- 
deGennes theory, which can be approximated with coupled time-independent 
gravitational Schrödinger equation. Further observation is of course recommended 
in order to refute or verify this proposition. 


1 Introduction 


In the past few years, there have been some hypotheses sug- 
gesting that quantization of planetary distance can be derived 
from a gravitational Schrédinger equation, such as Rubcic 
& Rubcic and also Nottale’s scale relativity method [1, 3]. 
Interestingly, the gravitational Bohr radius derived from this 
gravitational Schrödinger equation yields prediction of new 
type of astronomical observation in recent years, i.e. extra- 
solar planets, with unprecedented precision [2]. 

Furthermore, as we discuss in preceding paper [4], using 
similar assumption based on gravitational Bohr radius, one 
could predict new planetoids in the outer orbits of Pluto 
which are apparently in good agreement with recent observa- 
tional finding.. Therefore one could induce from this observ- 
ation that the gravitational Schrödinger equation (and gravi- 
tational Bohr radius) deserves further consideration. 

In the meantime, it is known that all present theories 
discussing gravitational Schrödinger equation only take its 
time-independent limit. Therefore it seems worth to find out 
the solution and implication of time-dependent gravitational 
Schrédinger equation (TDGSE). This is what we will discuss 
in the present paper. 

First we will find out numerical solution of time-inde- 
pendent gravitational Schrödinger equation which shall yield 
gravitational Bohr radius as expected [1, 2, 3]. Then we ex- 
tend our discussion to the problem of time-dependent grav- 
itational Schrédinger equation. 

In the subsequent section, we also discuss plausible ex- 
tension of this gravitational Schrédinger equation to include the 


effect of phion condensate via Gross-Pitaevskii equation, 
as described recently by Moffat [5]. Alternatively one can 
consider this phion condensate model from the viewpoint of 
Bogoliubov-deGennes theory, which can be approximated 
with coupled time-independent gravitational Schrédinger 
equation. To our knowledge this proposition of coupled time- 
independent gravitational Schrédinger equation has never 
been considered before elsewhere. 

Further observation is of course recommended in order 
to verify or refute the propositions outlined herein. 

All numerical computation was performed using Maple. 
Please note that in all conditions considered here, we use 
only gravitational Schrödinger equation as described in Rub- 
cic & Rubcic [3], therefore we neglect the scale relativistic 
effect for clarity. 


2 Numerical solution of time-independent gravitational 
Schrédinger equation and time-dependent gravita- 
tional Schrédinger equation 


First we write down the time-independent gravitational 
Schrödinger radial wave equation in accordance with Rubcic 
& Rubcic [3]: 


aR 2dR 8rm?EF' 








t R 
dr? rdr H2 T (1) 
24° GMM? p LEI) po 
7 H? r2 S 


When H, V, F' represents gravitational Planck constant, 
Newtonian potential, and the energy per unit mass of the 
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orbiting body, respectively, and [3]: 








H=n(an po), (2) 
mo 
GM 
V(r) =-——, 6) 
Hae. (4) 
m 


By assuming that R takes the form: 


R=e°" 


(5) 


and substituting it into equation (1), and using simplified 
terms only of equation (1), one gets: 


2ae 8nGMm? e72" 
r r H2 


After factoring this equation (6) and solving it by equat- 
ing the factor with zero, yields: 


2 (4nGMm? = H?a) 


—ar 
e,-ar 
e 


V=a 





(6) 





RR= ZH =0, (7) 
or 
RR =41GMm? — H7a=0, (8) 
and solving for a, one gets: 
4n?GMm? 


Gravitational Bohr radius is defined as inverse of this 
solution of a, then one finds (in accordance with Rubcic & 
Rubcic [3]): a 
~ 4n2GMm?’ 
and by substituting back equation (2) into (11), one gets [3]: 

2 2 
rı = (=) GM 3 
ac 
which is equivalent with Nottale’s result [1, 2], especially 
when we introduce the quantization number: rn = rın? [3]. 
For complete Maple session of these all steps, see Appen- 
dix 1. 

Solution of time-dependent gravitational Schrödinger 
equation is more or less similar with the above steps, except 
that we shall take into consideration the right hand side 


of Schrödinger equation and also assuming time dependent 
form of r: 


rı (10) 


(11) 


R=e or) | (12) 


Therefore the gravitational Schrédinger equation now 
reads: 








R  2dR. Snm?E! p 

dr? ` rdr ` H2 

2 4n°GMm* p é(é+1),_ ah ey 
r H2 r2 dt ’ 


or by using Leibniz chain rule, we can rewrite equation 
(15) as: 








dR dr(t) R 2dR | 8rm?E' p | 
dr(t) dt dr? ' rdr ` H? ! (14) 
2 2 
424 GMM’ p A 
r H? r2 


The remaining steps are similar with the aforementioned 
procedures for time-independent case, except that now one 
gets an additional term for RR: 


RR' = Ha Sr(t))r(e) — o’r(t)H? + (15) 
+8rGMm? — 2H?°a =0. 


At this point one shall assign a value for 4 r(t) term, 
because otherwise the equation cannot be solved. We choose 


4 r(t)=1 for simplicity, then one gets solution for (17): 








a=a,t7=7, m=m, H=H, G=G, M=M, 
a2:= 3 2 772 2 2 23/9 
t= RootOf(r(_Z)aH*—r(_Z)a H’ +8r GMm*—2aH*) 
{a=0, t=t, m=m, H=H, G=G, M=M, r=0}, 












































{a=0, r=7, t=t, m=m, H=H, M=M, G=0}, 

{r Tt, t=t, m=m, H=H, M=M,a=H,G atime} 
{a=a, H=0, r=n, t=t, m=m, M=M, G=G}, 
{a=0, r=n, t=t, m=m, H=H, G=G, M=0}, 
{a=0, r=7, t=t, H=H, G=G, M=M, m=0}, 

{a=a, H=0, r=7, t=t, m=m, G=G, M=0}, 

{a=a, H=0, r=7, t=t, G=G, M=M, m=0} 








Therefore one can conclude that there is time-dependent 
modification factor to conventional gravitational Bohr radius 
solution. For complete Maple session of these steps, see 
Appendix 2. 


3  Gross-Pitaevskii effect. Bogoliubov-deGennes appro- 
ximation and coupled time-independent gravitational 
Schrodinger equation 


At this point it seems worthwhile to take into consideration a 
proposition by Moffat, regarding modification of Newtonian 
acceleration law due to phion condensate medium, to include 
Yukawa type potential [5, 6]: 


GoM 











exp(—pgr 

a(r) = - +K ( 5 2 det ig (16) 
Therefore equation (1) can be rewritten to become: 
dR  2dR Sam? ED 
dr? ' rdr H? 
2 4r? (GM — K exp(—per)(1 + wgr))m? poy) 

r H? 
t+ pg 

r2 


or by assuming u = 2ko = Hor for the exponential term, eq- 
uation (17) can be rewritten as: 
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d’R p 2dR Brm?E' p 
dr? rdr | H2 i 
2 =) 2 
42 4n? (GM — Ke-?¥0(1+ por)) m = (18) 
r H2 
Et) p_ og 
72 


Then instead of equation (8), one gets: 
RR" =8rGMm?—2H’°a—8r’m? Ke” (1+u)=0. (19) 
Solving this equation will yield a modified gravitational 
Bohr radius which includes Yukawa effect: 
H? 
~ 4n?(GM — Ke-2#0)m? 





rı (20) 
and the modification factor can be expressed as ratio be- 
tween equation (20) and (11): 


GM 
a= 


= 21 
GM — Ke-?H0)’ en 


for complete Maple session of these steps, see Appendix 3. 

A careful reader may note that this “Yukawa potential 
effect” as shown in equation (21) could be used to explain 
the small discrepancy (around +8%) between the “observed 
distance” and the computed distance based on gravitational 
Bohr radius [4, 6a]. Nonetheless, in our opinion such an 
interpretation remains an open question, therefore it may be 
worth to explore further. 

There is, however, an alternative way to consider phion 
condensate medium, i.e. by introducing coupled Schrödinger 
equation, which is known as Bogoliubov-deGennes theory 
[7]. This method can be interpreted also as generalisation of 
assumption by Rubcic-Rubcic [3] of subquantum structure 
composed of positive-negative Planck mass. Therefore, 
taking this proposition seriously, then one comes to hypo- 
thesis that there shall be coupled Newtonian potential, in- 
stead of only equation (3). 

To simplify Bogoliubov-deGennes equation, we neglect 
the time-dependent case, therefore the wave equation can be 
written in matrix form [7, p. 4]: 





[A] [¥] =0, (22) 


where [A] is 2x2 matrix and [wv] is 2x1 matrix, respecti- 
vely, which can be represented as follows: 





81GMm7e%" ar 2ae %” 
rr nr 
[A] = (23) 
“ah. 20e AT 8rGMm? e72" 
ae 
Tr rh? 
and 
f(r 
[¥] = (r) . 24) 
g(r) 


Numerical solution of this matrix differential equation 
can be found in the same way with the previous methods, 
however we leave this problem as an exercise for the readers. 

It is clear here, however, that Bogoliubov-deGennes ap- 
proximation of gravitational Schrödinger equation, taking 
into consideration phion condensate medium will yield non- 
linear effect, because it requires solution of matrix differen- 
tial equation* (22) rather than standard ODE in conventional 
Schrödinger equation. This perhaps may explain complicated 
structure beyond Jovian Planets, such as Kuiper Belt, inner 
and outer Oort Cloud etc. which of course these structure 
cannot be predicted by simple gravitational Schrödinger eq- 
uation [1, 2, 3]. In turn, from the solution of (22) one could 
expect that there are multitude of celestial objects not found 
yet in the Oort Cloud. 

Further observation is also recommended in order to 
verify and explore further this proposition. 


4 Concluding remarks 


In the present paper, a numerical solution of time-dependent 
gravitational Schrödinger equation is presented, apparently 
for the first time. This numerical solution leads to gravita- 
tional Bohr-radius, as expected. 

In the subsequent section, we also discuss plausible ex- 
tension of this gravitational Schrödinger equation to include 
the effect of phion condensate via Gross-Pitaevskii equation, 
as described recently by Moffat. Alternatively one can con- 
sider this condensate from the viewpoint of Bogoliubov- 
deGennes theory, which can be approximated with coupled 
time-independent gravitational Schrödinger equation. 

It is recommended to conduct further observation in order 
to verify and also to explore various implications of our pro- 
positions as described herein. 
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Appendix 1 Time-independent gravitational Schrödinger equation 


> restart; 
> with (linalg); 
> R: = exp(—(alpha*r)); 
R = e T 


> DIR:=diff(R,r); D2R:=diff(D1R,r); 


D1R:= -ae ®" 
D2R := —a*e%" 


> SCHEQ1:=D2R+D1R*2/r+8*pi"2*m*E*R/h`2+8*pi^2*G*M*m^2*R/(r*h^2)— 
1*(1+1)*R/r^2=0; 
> XX1:=factor (SCHEQ1); 
> #Using simplified terms only from equation (A*8, of Rubcic & Rubcic, 1998) 
> ODESCHEQ:=D2R+D1R*2/r+8*pi’2*G*M*m’2*R/(r*h’2)=0; 
Qae*" 8n2GMme%" 


ODESCHEQ :=a°e *" =0 
r rH? 





> XX2:=factor (SCHEQ2); 
e72" (a?rH? — 2H°a + 8n*GMm?) 
PH? 


XXQ2:= =0 





> RR:= solve (XX2, r); 
2(4n?GMm? — H?a) 


RR: a2 2 








> #Then solving for RR=0, yields: 
> SCHEQ3:= 4*pi^2*G*M*m^2—h^2*alpha=0; 


SCHEQ3 := 4r? GMm? — H?a=0 


> a:= solve (SCHEQ3, alpha); 
E 4n?GMm? 
oe 
> #Gravitational Bohr radius is defined as inverse of alpha: 
> gravBohrradius:= 1/a; 
H? 
4n2GMm? 


TgravBohr := 


Appendix 2 Time-dependent gravitational Schrodinger equation 


> #Solution of gravitational Schrodinger equation (Rubcic, Fizika 1998); 
> restart; 
> #with time evolution (Hagendorn’s paper); 
> S:=r (t); R:=exp(—(alpha*S)); R1:=exp (—(alpha*r)); 
S := r(t) 
R — e~ r 


> D4R:= diff(S,t); DIR:=—alpha*exp (—(alpha*S)); D2R:=—alpha^2* 
exp (—(alpha*S)); D5R:=D1R*D4R; 


d 
D4R := — r(t) 
dt 


DIR := -a e7% "T®) 


D2R := —a? e7% "®) 


d 
DIR := -ae~*") = r(t) 
dt 


> #Using simplified terms only from equation (A*8) 
> SCHEQ3:= —h*D5R+D2R+D1R*2/S+8*pi'2*G*M*m’2*R/(S*h’2); 
> XX2:=factor (SCHEQ3); 


eT ort) (Hat r(t)—a?r(t)H?-2H? a+ 8n?GMm?) 


r (t)H? ze 


XXQ = 





> #From standard solution of gravitational Schrodinger equation, we know (Rubcic, 
Fizika 1998): 
> SCHEQ4:= 4*pi*2*G*M*m’2—h’2*alpha; 


SCHEQ4 := 4n°GMm? — H?a 


> #Therefore time-dependent solution of Schrodinger equation may introduce new 
term to this gravitational Bohr radius. 
> SCHEQS:=(XX2*(S*h°2)/(exp (—(alpha*S))))—2*SCHEQ4; 
dr(t 

ODESCHEQ5 := H’ a ar) 5) — a? r(t)H? 
> #Then we shall assume for simplicity by assigning value to d[r(t)]/dt: 
> D4R:=1; 
> Then we can solve again SCHEQS similar to solution of SCHEQ4: 
> a2:=solve((h*3*alpha*(D4R)*S —alpha’2*S*h’2)+2*SCHEQ4); 








a=a,T=7,m=m, H=H, G=G, M=M, 
a2 :=\ t = RootOf(r(_Z)aH?—r(_Z)a?H?+82?GMm?—2aH?) } 















































{a=0, t=t, m=m, H=H, G=G, M=M, r=0}, 

{a=0, r=7, t=t, m=m, H=H, M=M, G=0}, 

{r T, t=t, m=m, H=H, M=M, a=H, G= 5 ie 
{a=a, H=0, r=7, t=t, m=m, M=M, G=G}, 

{a=0, 7 , t=t, m=m, H=H, G=G, M=0}, 

{a=0, r=7, t=t, H=H, G=G, M=M, m=0}, 

{a=a, H=0, r=7, t=t, m=m, G=G, M=0}, 

{a=a, H=0, R= Tt t, G=G, M=M, m=0} 





> #Therefore one could expect that there is time-dependent change of gravitational 
Bohr radius. 


Appendix 3 Time-independent gravitational Schrödinger equation 
with Yukawa potential [5] 


> #Extension of gravitational Schrodinger equation (Rubcic, Fizika 1998); 
> restart; 
> #departure from Newton potential; 
> R:=exp (—(alpha*r)); 
R:=e “" 


> DIR:= diff(R,r); D2R:= diff(D1IR,n); 


DIR :=-ae *" 


2,-ar 


D2R := -—a“e 
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> SCHEQ2:=D2R+D1R*2/r-+8*pi’2*(G*M—K*exp (—2*mu)*(1-+-mu*r))*m’2*R/ 





=0 





(*h’2)=0; 
2 — cr 
ODESCHEQ := ae *" — “2° 
r 
; 87r? (GM — Ke™?H(1 + ur)jm? e72" 
i 


rH? 
> XX2:=factor(SCHEQ2); 
> RRI1:=solve(XX2,r); 
2(—H?a + 4n?GMm? — 422m? Ke~?+) 
—a? H? + 812m? Ke-72. 


RR1: 








> #from standard gravitational Schrodinger equation we know: 
> SCHEQ3:=4*pi'2*G*M*m‘*2—h’2*alpha=0; 
> a:=solve(SCHEQ3, alpha); 
> #Gravitational Bohr radius is defined as inverse of alpha: 
> gravBohrradius:= 1/a; 

H? 
4n?GMm? 


TgravBohr >= 


> #Therefore we conclude that the new terms of RR shall yield new terms (YY) into 
this gravitational Bohr radius: 

> PI:= (RR*(alpha’2*h*2)—(—8*pi'2*G*M*m‘°2+2*h*2*alpha)); 

> #This new term induced by pion condensation via Gross-Pitaevskii equation may 
be observed in the form of long-range potential effect. (see Moffat J., arXiv: astro- 
ph/0602607, 2006; also Smarandache F. and Christianto V. Progress in Physics, v. 2, 
2006, & v. 1, 2007, www.ptep-online.com) 

> #We can also solve directly: 

> SCHEQS:=RR*(alpha*2*h’2)/2; 


SCHEQS := a? H?(—H?a + 4n?GMm? — 4n?m? Ke~?) 
l —a? H? + 8r?m? Ke-2# 





> al:=solve(SCHEQS, alpha); 





4n?m?(GM — Ke~?) 
al:=0,0, 
H2 
> #Then one finds modified gravitational Bohr radius in the form: 
> modifgravBohrradius:=1/(4*pi"2*(G*M—K *exp (—2*mu))*m^2/h^2); 
H? 
4n?2m?(GM — Ke-24) 





Tmodified.gravBohr ‘= 
> #This modification can be expressed in chi-factor: 
> chi:=modifgravBohrradius/gravBohrradius; 


GM 


X= GM — Ke? 
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